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Abstract. The concept of soliton complex in a nonlinear dispersive medium is pro- 
►^ ] posed. It is shown that strongly interacting identical topological solitons in the medium 

O^ ' can form bound soliton complexes which move without radiation. This phenomenon 

f^ ' is considered to be universal and applicable to various physical systems. The soliton 

f^ , complex and its "excited" states are described analytically and numerically as solutions 

O^ ' of nonlinear dispersive equations with the fourth and higher order spatial or mixed 

^^ ' derivatives. The dispersive sine-Gordon, double and triple sine-Gordon, and piecewise- 

0> , linear models are studied in detail. Mechanisms and conditions of the formation of 

soliton complexes, and peculiarities of their stationary dynamics are investigated. A 
phenomenological approach to the description of the complexes and the classification 
of all the possible complex states are proposed. Some examples of physical systems, 
where the phenomenon can be experimentally observed, are briefly discussed. 
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I. INTRODUCTION 

Properties of wave excitations in condensed matter are strongly influenced 
by spatial dispersion. In solids dispersion originates from the discreteness of real 
crystals. The nonlinear dynamics of lattice models exhibits thus many specific 
phenomena known as the discreteness effects [1 — 7]. In the long wave limit, 
when the systems are considered as continuous, these effects have to disappear. 
However some of them leave a trace in the continuum limit, and hence can be 
picked as a class of universal dispersive effects. 



Naturally the same situation is observed in macroscopic discrete systems, 
e.g., arrays of Joscphson junctions and nonlinear transmission lines, systems 
with nonlocal interactions, and others [8 — 15]. 

The universal effects can also manifest themselves in strongly dispersive 
media, in which the dispersion is not a consequence of translational symmetry 
of underlying structures. Examples of such systems are plasmas, fluids, optical 
and dissipative-dispersive systems [16 — 23]. 

One of these effects, occuring in the dynamics of nonlinear excitations of a 
dispersive medium is the formation of bound states of solitons [9, 24] . The role of 
dispersion as a factor influencing the interaction between well-separated solitons 
was discussed in a great number of works (see, e.g., [9, 18, 22 — 24]). As a result, 
two mechanisms of a formation of soliton bound states were found out. In the 
first case interactions of oscillating soliton tails in dispersive media lead to the 
formation of bunches of solitons, consisting of two or more well-distingushable 
humps [18,22]. In the second case solitons coexisting with resonant radiation 
can form bound states with purely solitonic asymptotics due to some kind of the 
radiation interference effect [25 — 27]. In these theories the contribution of the 
dispersion to the soliton interaction has to be considered as a weak perturbation. 
As a result the energy of the formed multisoliton structure differs slightly from 
the energy of a corresponding set of the non-interacting solitons. 

In the case of topological solitons this phenomenon of soliton bunching was 
observed through numerical simulations beginning with the work [8], where, 
in fact, the weakly discrete sine-Gordon model was studied. Then there were 
attempts [28, 29] at explaining the effect basing on the use of the soliton per- 
turbation theory applied to the continuous analogue of the system. Authors of 
[28] were the first who pointed out that any dispersive terms should be added to 
the usual sine-Gordon model to obtain the multisoliton steady solutions. Then 
Peyrard and Kruskal observed by means of a numerical simulation the almost 
radiationless motion of the 47r-soliton in the highly discrete sine-Gordon sys- 
tem [2]. This effect occurs not only in the discrete model [2,8, 15] but also in 
the continuous dispersive sine-Gordon equation, and a corresponding analytical 
solution for the 47r-soliton complex can be found exactly [30, 31]. 

In general, topological multisolitons exist as discrete sets of solitonic con- 
figurations with internal structures. They are well-studied in systems with 
nonlocal interactions [11,12,14]. At last, bound solitons can also be realized 
in systems of anharmonically interacting particles. Such models are described 
by the Boussinesq-type equations with the high spatial derivatives [22,32,33], 
and properties of the bound solitons in these dispersive systems were studied in 
detail in both cases of topological and non-topological solitons [22, 32 — 37]. 

In the present paper we concentrate just on the case of strong interaction 
between topological solitons. This occurs when identical solitons are closely 
placed. As a result the repulsive potential of the solitons has to grow rapidly 
with decreasing the distance between them. This causes the applicability of 
the perturbation theory, basing on the " one-soliton" approximation to become 



invalid. Then one has to take into account the dispersion, as an additional 
source of strong interaction, and this can change the character of the interaction 
between solitons in general. 

The aim of this paper is to show that solitons in a strongly dispersive medium 
possess an internal structure and their interaction depends on intrinsic proper- 
ties such as flexibility. Due to this dependence the potential energy turns out to 
be a non-monotonic function of the distance between solitons. As a result identi- 
cal solitons can attract each other and form a bound-soliton complex which can 
move without any radiation in strongly dispersive media [30,31]. Thus we call 
the bound soliton states with the zero and small distances between composite 
solitons the soliton complex and its "excited" states respectively. 

We present a number of dispersive models bearing such topological soliton 
complexes. The models are described by nonlinear equations with fourth and 
higher spatial or mixed derivatives. Solutions of the relevant equations can 
be obtained numerically and analytically. We found exact analytical solutions 
for two variants of the dispersive sine-Gordon (dSG) and double sine-Gordon 
equations (dDSG). We aslo show that the complex, consisting of three solitons, 
can be described explicitly in the dispersive equations with the additional sixth 
derivative in the cases of the sine-Gordon and triple sine-Gordon models. 

We propose the classification of the "excited" states of soliton complexes, 
constructing them explicitly in the framework of the double piccewise-linear 
dispersive model. The two-soliton ansatz approximation is used to establish 
analytically the existence condition for the soliton complex in the dSG and 
dDSG equations. We find numerical solutions of these equations for the two- 
soliton complex and its "excited" states, and their dependences of energies and 
velocities on the model parameters. As a result we formulate the concept of the 
soliton complex and classify it as a specific bound state of strongly interacting 
identical solitons in a dispersive medium. 

2. DISPERSIVE MODELS WITH SOLITON COMPLEXES 

As a first example of a dispersive model we mention the discrete sine-Gordon 
system which is described by the equation [2, 4, 38]: 

— ^ -f- 2Un - Un-l - Un+1 + -^ Sm(u„) = 0, (1) 

where u„ is, e.g., the displacement of atom n and d is the discreteness parameter. 
A stationary motion of a single 27r-soliton is impossible in this dispersive system 
because of a strong radiation emitted by the moving 27r-soliton [4] . At the same 
time numerical simulations [2, 8, 15] showed the almost radiationless motion of 
the 47r-soliton and other 2?Ti7r-soliton complexes. Authors of [2] tried to explain 
the formation of the soliton complex of two identical 27r-solitons by exploiting 
the fact of the presence of the Peierls potential in the lattice under consideration. 



However in works [30, 31] it was found that the radiationless motion of the 
complex can be described explicitly in the framework of the dispersive sine- 
Gordon equation (IdSG) with a fourth spatial derivative: 

utt - u^x - (iuxxxx + sin(u) = 0. (2) 

Eq. (2) is obtained as the long wave limit of Eq. (1) by substituting xd for n, 
td for r, and the second difference by the series: 

Uji—i -\- Uji^i ZUn ~ Uxx ~r ptlxxxx ~r ■■■ \0) 

The relation between the dispersive factor (3 and the discreteness parameter d 
is /3 = l/(12(i^). The exact solution of Eq. (2), corresponding to the soliton 
complex, has the following form: 

M47r = 8arctan{exp(W-— ^^=)}, (4) 




yo = ±\i-\/^. (5) 



The velocity of the complex is not an arbitrary constant in the solution (4), but 
it is a function of the dispersive parameter (3. As a function of the discreteness 
parameter, it equals Vo{d) — ±-\/l — (l/Sd) and differs less than five percents 
from numerical result of Peyrard and Kruskal, as it follows from a comparison of 
Fig. 1 (solid line) and Fig. 12 of [2]. This fact provides the starting point of our 
investigation because it shows that taking into account a high-order dispersion 
in the continuum model described by Eq. (2), leads to the same phenomenon 
as that in the discrete model of Eq. (1). 

From the other hand it is clear that the continuum model of Eq. (2) describes 
properly stationary moving nonlinear excitations in the original discrete model 
of Eq. (1) only in the limit (3^1. In this case using Eq. (4) one can write the 
discrete 47r-soliton in the first approximation as 

p}, 
Uji — 8arctan{expfW— (n — Vqt)]}. 
V d 

For small /3 (large d) the solution is a smoothly-varying function of number n, 
and its effective width is proportional to yd. This provides a validity of the 
continuum approximation. It is remarkable and surprising that the analytical 
expression Vo('^) for a velocity holds good in the highly discrete system, i.e., far 
from the continuum limit. It would be also noted that the 47r-soliton appears to 
be localized more strongly than the 27r-kink of the usual sine-Gordon equation, 
since the kink width is proportional d. 

It is evident [31] that solutions similar to Eq. (4) are available for the second 
or "regularized" dispersive sine-Gordon equation (2dSG) which has the fourth 



spatio-temporal mixed derivative instead the fourth spatial derivative of Eq. 

(2): 

utt - u^x - (iuttxx + sin(u) == 0. (6) 

The form of the soliton complex solution of the equation is the same as for the 
dSG equation, however the velocity dependence on the parameter (5 differs from 

Eq. (5). 



K(^)=±(Vi+f-vf)- (^) 

As a function of the parameter d the velocity is also shown in Fig. 1 (dash line). 
The soliton complex with the internal structure was virtually identified by 
Peyrard et al [39] in simulations of the continuous modified sine-Gordon model. 
In general, the taking into account the spatial dispersion or nonlocal interac- 
tions leads to a possibility of the complex formation [40] . In fact, the presence of 
freely moving soliton complexes and their " excited states" was examined numer- 
ically in the continuous nonlocal sine-Gordon models in works [11, 12, 14]. Such 
models are described by integro-differcntial equations which can be reduced to 
systems of two local equations in the case of exponentially-decaying kernel. As 
an example, we point out the equation [11]: 

oo 

d f 

Utt + sm{u) ^ — / dx'G{x — x')uxi{x' ^t). (8) 

— oo 

In the case of the kernel G(a;, A) = (1/2A) exp(— |a;| /A) this equation is trans- 
formed into the set of equations 

\^Wxx-W^-Ux, Mtt +Sin(u) = Wa;, 

which possess soliton-complex solutions. 

All the above facts collected together prompt that the formation of the soli- 
ton complexes is a universal property of the strongly dispersive media. There- 
fore, one can try to find solutions for the soliton complexes in more general 
systems than the dispersive sine-Gordon models. 

The present paper is mainly devoted to studying the soliton complex in the 
dDSG equations because of two reasons. First, the usual DSG equations describe 
a large variety of physical systems: ferro- and antifcrromagncts, magneto-elastic 
systems, superfiuid "^He and others [41 — 43]. Secondly, in the usual DSG equa- 
tion 27r-kinks form the 47r-kink, the wobbler, due to the action of an external 
field. In the dispersive DSG equations both factors, dispersion and external 
field, cause the soliton coupling, and it is interesting to investigate their mu- 
tual influence on the binding process. At last the dDSG equations contain the 
sine-Gordon equations as the limit cases. 



Thus we deal with the dispersive double sine-Gordon equation (IdDSG): 

utt - Uxx - (3uxxxx + sm(u) + 2hsm{-) = 0, (9) 

and its regularized variant (2dDSG): 

u 
Utt - Uxx - ffuttxx + sin(u) + 2/isin(-) = 0, (10) 

where, for example, in magnetic applications 4>{x,t) = ^u{x,t) denotes the 
azimuth angle of the magnetization vector in the easy-plane ferromagnet, and 
/i is a magnetic field applied along the easy plane. When h = the equations 
revert to the dispersive and regularized sine-Gordon equations. 
Eq. (9) can be derived from the Lagrangian: 

oo 

/I U 

-K -ul + (iul, - 2(1 - cos{u)) - 8/1(1 - cos{-))}dx. (11) 



The Hamiltonian of the dispersive double sine-Gordon system is given by: 



H 



/■ 1 u 

J ^Wt + ul - Pul^ + 2(1 - cos(m)) + 8/i(l - cos(-))}cia:. (12) 



Corresponding expressions for the regularized Eq. (10) arc obtained by substi- 
tuting Mj^, for u"^^ in Eqs. (11) and (12). 

The difference between the dDSG equations results in two types of spectra 
of linear excitations. Namely, Eq. (9) and Eq. (10) have the dispersion relations 
u}{k) = ^1 + h + k^ - (3k* and uj{k) = .^{1 + h + k^)/{l + pk^) , respectively. 
For the first spectrum there exists formally the critical wave number fco at which 
'^(^o) = . This means that equilibrium state u = is unstable with respect to 
the short wave length perturbations. Recalling the spectrum of linear excitations 
of the discrete system one realizes the artificial origin of this instability. The 
regularized equation is introduced in order to avoid this instability [33, 44] as 
easily seen from its spectrum. Moreover, such a spectrum shares the main 
features of Eq. (8), the sine-Gordon model with the nonlocal interaction. 

Some important properties of the dSG and dDSG equations can be repro- 
duced by the following model equation which we call as the dispersive double 
piecewise-linear equation: 

Utt - Uxx - (iUxxxx + .f{u) = 0, (13) 

where force f{u) is a periodic function of the period 4tt. On the interval [—2tt, 27r] 
it is given by 



(1 + /i)(u + 27r) -2tt < u < -uq 
f{u) -<((!- h)u \u\ < uo (14) 

(1 + h){u - 2tt) uq < u < 27r. 

The double quadric potential of the model is presented in Fig. 2. The point, 
corresponding to the energy niaxinium, is chosen as uq = 2 arccos(/i) to model 
the behaviour of the double sine-Gordon model. As h — and h = 1 the 
equation degenerates into the analogue of the dSG equation with the periods 
2tt and An, respectively. We show further that this model allows to construct 
explicitly the expressions for the soliton complexes and their "excited" states. 

Other dispersive equations combining the properties of the Klein-Gordon and 
Boussinesq-type equations are known [34 — 36] , in which topological solitons 
are shown analytically and numerically to possess an internal structure, and, 
therefore, they are good candidates for bearing soliton complexes. 



EXACT SOLUTIONS FOR SOLITON COMPLEXES AND 
THEIR DYNAMICAL PROPERTIES 



The exact moving solution of Eq. (4) describes the soliton motion with- 
out radiation in a dispersive medium. This radiationless dynamics of a soliton 
complex in the dSG equation was discussed in detail [30,31]. In this section 
we present exact solutions for the dispersive double sine-Gordon equations and 
discuss the peculiarities of their dynamics. For the concrete definition we imply 
the magnetic applications of the equations, where solitons are domain walls and 
parameter /i is a magnetic field. Therefore we use further this terminology. 

It is known that applying a magnetic field to a ferromagnet leads to coupling 
180° domain walls into 360° domain wall. In terms of the usual double sine- 
Gordon equation (/3 = 0) this means the existence of a wobbler solution [41, 42] 
or the 47r-kink. Such a static solution does not exist for Eq. (9) if /3 7^ but it 
holds for Eq. (10): 

■Uto(x) — 4arctan{exp((7tuX — Rw)} + 4arctan{exp((7tuX + Rw)} (15) 

where q^ = \/(l -I- h) and sinh(i?tu) — l/^/h . 

In this work we are interested in a stationary soliton motion, i.e. consider 
the solutions of the form u{x, t) = u{x — Vt). Then both the equations (9) and 
(10) are reduced to the ordinary differential equation: 

u 
u^^ + auzzzz - sin(M) - 2ft,sin(-) = 0. (16) 

Here z — [x — Vt)/^/l — V'^, and parameter a equals to a'-^-' (a'-^-') for the 
IdDSG (2dDSG) equation, where 



(1) P (2) Py 



-''' - TT^TT^^ -''' - TT^^T^.- ^ 



When /i = 0, Eq. (16) is reduced to the dSG case which was analysed in [30, 31]. 
Some results of the numerical integration of the reduced equation, as the specific 
limit of the nonlocal model (8), were discussed in [11]. 

When /i 7^ we are able to find again the exact solution of Eq. (16): 



/2 x—Vt 

UA-n = 8 arctan{exp(W - + h )}. (18) 

Thus the IdDSG and 2dDSG equations have complex solutions in identical form 
but with different expressions for the velocities 



yi(/3,/i)-±Wl-(2 + 3/i)\/^. (19) 




y2(/3, M = ±(^1 + f (1 + \hY - y f (1 + \h)). (20) 

These solutions describe complexes consisting of two strongly bound 27r-solitons. 
They differ from the wobbler Eq. (15) by the effective widths, the zero distance 
between solitons, and the ability to move in the dispersive medium. Since now 
velocities are functions of parameter /i, we can change them from the maximum 
values corresponding to the dSG limit (Eq. (5) and Eq. (7)) to zero. In Eq. (9) 

it occurs at the finite critical value her = \/ 573 ~ f • Fo'^ example, in ferromagnets 
we can control the velocity of a motion of the domain wall complex through the 
magnetic field. 

The next peculiarity of the complex dynamics is revealed, if we evaluate 
the energies using the Hamiltonian expressions (see Eq. (12)). For the IdDSG 
complex we obtain: 

ii;i^32[(3/3)-^-^(3/3)^]. (21) 

The energy of the soliton complex turns out to be independent of the parameter 
ft., and, therefore, of the velocity at the fixed /?! This means that, at least, by 
abiabatic variation of the magnetic field, we can vary the form and speed of the 
complex, conserving its energy. This remarkable property could be used in the 
energy transfer applications in physical systems described by the IdDSG equa- 
tion. However, the property is very much sensitive to spectrum characteristics 
of the dispersive medim. The regularized Eq. (12) has the velocity and field 
dependent energy: 

E^ = 32[(3/3y2')-^ - ^(3/31^2')^], (22) 



where V2{P,h) is given by the expression (20). More details of comparative 
analysis of the energy dependences of the complexes are presented in the section 
6. 

The next natural question is the existence of an exact solution for soliton 
complexes consisting of more than two solitons. Numerical integrations of the 
discrete sine-Gordon equation Eq. (1) and the nonlocal model Eq. (8) re- 
veal such solutions. We are able to find the analytical solution describing the 
three-soliton complex in the dispersive sine-Gordon equation with sixth spatial 
derivative [45]: 

utt - Uxx - Puxxxx - lUxxxxxx + sin(M) = 0. (23) 

The equation can be derived from a discrete model if we take into account the 
higher-order term in an expansion of the type of Eq. (3) . For the special choice 
of the parameters 7 = Jj/3 2 , Eq. (23) has the following exact solution: 




Mstt ^ 12arctan{cxp(^/— )}, (24) 

45 ^1 - v;2 



where the velocity takes the form: 



K^±f-|V^. (25) 

The exact solution of the form of Eq. (24) exists also in the dispersive triple 
sine-Gordon equation [45]: 

u 2 

Utt - Uxx - (iuxxxx - JUxxxxxx + sin(u) + hi sin(-) + /12 sin(-u) = 0, (26) 

which is the generalization of Eq. (23), and where parameters hi and /12 are 
arbitrary constants. Under the condition 7 = ^/3^(1 + ^^2) the thrce-soliton 
complex is described by the expression: 

u^TT — 12arctan{cxp(g3(x — Vst)}. (27) 

The relations between solution parameters are the following: 

'/a = (^)^/^ (28) 




V, ^±dl-{- + hi- -h2)J-^—. (29) 



In previous formulas for the presentation of exact solutions we have used the 
Lorentz-invariant-like expressions to keep the analogy with those of the Lorentz- 
invariant equations with (3 = 0. However the solutions can be written also in 



the form like Eq. (27) which is evidently simpler. In any case at a first glance 
the solutions for the soliton complexes and conditions of their existence look 
like exotic. This raises, at least, three questions. The first of these: how are 
the solutions sensitive to the variation of the equation parameters? The second, 
what kind of other solutions exists in the strongly dispersive equations? And 
the third, what is the mechanism underlying the creation of these bound states? 
Pondering these questions provides the substance of the next sections. 



4. PHENOMENOLOGY OF SOLITON COMPLEX 
A. Collective coordinate approach to soliton-complex formation 

An analytical approach to the description of the soliton-complex formation in 
the dSG equation was proposed in [30, 31] . It is based on the use of the collective 
variable ansatz which is constructed by taking into account the translational and 
internal degrees of freedom of a soliton as well as interactions between solitons 
and solitons with radiation: 

u{x,t) =M('')(x,t;/,X,i?)+u('')(x,t). (30) 

Here m^^^ is a solitonic part and u^''''{x, t) is a part of the solution corresponding 
to radiation. It turns out that the condition of the complex formation of the 
closely sited solitons can be found from the energy expression of the pair of 
strongly interacting solitons without taking into account the radiation [31]. Now 
we use this approximation for the description of the dispersive double sine- 
Gordon system. So we suppose that the complex dynamics can be considered 
in the framework of the soliton ansatz: 

X — X X — X 
y(s) ^ 4 arctan{exp( — R)} + 4 arctan{exp( — hi?)}, (31) 

which is prompted by the forms of the wobbler Eq. (15) and the exact solution in 
Eq. (18). Here I = l{t), X = X{t) and R = R{t) are functions of time. Functions 
l{t) and X{t) describe the changing of the effective width of solitons and their 
translational motion, respectively. The function R = R{t) corresponds to the 
changing separation between solitons, which is defined obviously as L = 21R. 

Inserting the ansatz into Eq. (11) and Eq. (12) we find the effective La- 
grangian and Hamiltonian of two interacting solitons in strongly dispersive me- 
dia. 

Leff^T-U, H,ff = T + U, (32) 

where T and U , the kinetic and potential energies, are given as 
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l;.2ron2 I /'r?2 , TT 



T - f{^l/[2R^ + (i?2 + \)h{R)] + {lRttSiiih{R)yi2{R) + 2RRtllt ^.^^. 

+Xt^h{R)}, 



(34) 



and 

U ^ f {/i(i?) + PhiR) + 2hPloiR) 
^f[l-tanh^(i.)(l-/.(i?)(l + -^))]}, 

where the foUowing notations have been introduced (see also Appendix A): 

/o(i?) = 2i?coth(i?), 

^i(^) = ^+sinhf2i?)' 
/.(i?)=coth^(i?)(l-^^^). 

Details of the calculation of explicit expressions for the kinetic and potential 
energies are presented in Appendix A. The expressions look rather complicated, 
but they contain a rich information about two-soliton dynamics and interactions. 
In particular, in the limit cases they describe the exact two-soliton solution of 
the integrable sine-Gordon equation and the small wobbler oscillations in the 
usual double sine- Gordon equation (see Appendix B and C). 

Since we are interested in the stationary soliton complexes let us investigate 
at first the potential energy of two interacting solitons. The energy is drawn 
in Figs. 3(a)-3(d) as a function of parameters I and R for two different values 
of the dispersive parameter /?. For the sake of symmetry we show the energy 
dependence on both positive and negative values of R, because the ansatz (31) 
and the energy are even functions of R. 

Let us examine the character of the dependence of the potential energy on 
the separation between solitons at the various /. In Figs. 3(a) and 3(b) we take 
P — Pi — ^. One can see that in the whole area shown the closely placed solitons 
repulse each other and the energy has a maximum at i? = 0. There exists the 
critical effective length l^r at which two local minima appear for the first time 
at i? := oo, and then, when / increases, they quickly diminish their separation. 
The minima correspond to the equilibrium positions of two solitons in which 
their mutual repulsion is balanced by the opposite action of the dispersive part 
of the interaction and the magnetic field. We show two cases, when h — and 
ft, = 0.1 with a view to compare results for the sine-Gordon and double sine- 
Gordon systems. Fig. 3(b) demonstrates the strong attractive contribution of 
the magnetic field to the soliton interaction. 

In Figs. 3(c) and 3(d) the parameter (3 is chosen as /3o = f- One can find that 
in this case the critical value Icr corresponds the point where the local maximum 
in the i?-dependence is changed by a local minimum. This obviously means that 
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in this point the repulsion between solitons is changed into the attraction. This 
fact can be proved strictly and the critical values can determined exactly because 
at small R the consideration can be performed analytically (see Appendix A for 
details). Indeed, from the Largangian expression Eqs. (87-89) valid for small R 
one derives the following Langrange equations: 

^llu ^;? - 1 + Xf + (i + h)P + ^ = 0, (35) 

—— ^ p = const, (36) 

Xf-l + i^ + h)f+^-^^0. (37) 

The stationary values of Iq and Xt = Vq are easily found. They coincide with 
parameters of the exact solution for the soliton complex, Eqs. (18) and (19), 
and give the critical values for the parameters at which the soliton attraction 
arises. It is easy to see, when Eq. (37) is satisfied, that the contribution of 
terms of O(R^) to the "stationary" Lagrangian part Eq. (89) equals zero, i.e. 
at this point it changes sign. That is, for example, for Vb = and h = 0, 
one finds from Eqs. (19) the following critical parameter values: /3o = | and 
'o(/3o) — Icr — ^73/2 w 1.225. In fact, at this value Icr the potential energy 
exhibits the flat plot at small R as it is seen in Fig. 3(c). Hence the possibility 
of solitons to change their effective length in the dispersive medium leads to 
changing the very character of interaction between them. Thus the exact com- 
plex solutions correspond to the bound states of solitons coupled in reason of 
their own attraction. 

Stability of the complex with respect to changing its parameters, velocity 
and effective width, as well as a possibility of its decay can be tested by the 
numerical simulation of the equations. This work is in progress. Results of 
[2, 11, 15] give information about the stability of the complexes in the discrete 
and nonlocal sine-Gordon model. The complexes manifest themselves as attrac- 
tors in the soliton dynamics of these dispersive media. Two solitons, moving 
with velocities larger than the critical value radiate energy until the velocity 
reaches its stationary value. Then the formed complex moves radiationlessly. 
At a velocity smaller than the critical value solitons repulse each other, and the 
complex decays. As a result the composite solitons travel from the center to 
their new equilibrium positions. 

This picture is consistent with conclusions following from our direct energy 
analysis. Further consideration shows that, for the first time, equilibrium local 
minima at _R = oo appear for the following values of the parameter I: 



m = \j{[\+p{i + h)Y'^ + \}i{i + h). 



(38) 
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It is easy to see that the parameter lo of the exact complex solution also be- 
longs to this dependence. Therefore at the same moment , when the attraction 
between solitons arises, the equilibrium local minima disappear and vice versa. 
The presence of the local minima points out on a possibility of the existence 
of another attractor of the wobbler-like type in soliton dynamics. To study 
the question about the existence of other stationary moving bound states of 
identical solitons, besides the exact soliton complex, we recur to the ordinary 
differential equation (16). 

B. Excited states of the soliton complex 

The equation (16) for stationary states has the asymptotics u(z) = Acxp{qz) 
where q obeys the equation 

aq^ + q^ - I - h = 0. (39) 

Soliton solutions are characterized by the vanishing asymptotics, hence the cor- 
responding q — K is supposed to be a real parameter: 



n=^{[^+a{l + h)]^/^-^}/a. (40) 

The exact complex solution is realized when a — 3/4. 

However Eq. (16) is of the fourth order and has another asymptotics which 
is oscillating one. In this case q — ik, where 



k^^{[\ + a{l + h)V/^ + ^}/a. (41) 

In general, a solution of Eq. (16) includes the solitonic and oscillating parts. 
This is simply interpreted, because a single moving 27r-soliton in the dispersive 
system usually radiates energy and generates the continuous waves. Thus in 
the case of the moving soliton we deal with the self-modulated medium. It is 
interesting to note that it occurs independently of the stability property of the 
equilibrium state. In fact, both equations (9) and (10) are reduced to Eq. (16), 
in spite of different dispersion relations. 

When two solitons are present in the system, they create the radiation back- 
ground, and one soliton moves upon the undulations generated by the other 
soliton. The situation turns out to be similar to that in the original discrete 
system [2], where solitons travel upon the Peierls potential. The existence of the 
purely solitonic stationary excitations under such conditions implies that some 
interference effect takes place, which leads to cancelling the radiation far from 
soliton complex. Available theories [25-27] of a soliton binding by radiation 
field suggest a large separation between composite solitons and a small influ- 
ence of the dispersion on soliton interactions. It is not correct for closely-sited 
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solitons forming the soliton complex. As a result an alternative approach to the 
description of these bound states is required. 

To find all possible forms of the soliton-complex solutions one has to solve 
the nonlinear eigenvalue problem, such as Eq. (16), which is hardly feasible by 
analytical tools. There are known some attempts to solve similar problems by 
variational methods [26,46]. We propose another approach to a searching for 
soliton-complex solutions and demonstrate it by applying to Eq. (16). 

So we seek solutions of Eq. (16), imposing the vanishing boundary condi- 
tions. Let us reformulate this problem as a self-consistent eigenvalue problem 
for the linear equation: 



dz 
where the potential well is 



[-"73 + U{z)]u,, = 0, (42) 



U{z) = '^^(^) + ^hMu/2) ^ 

Uzz 

We know that Eqs. (42) and (43) have, at least, the solution corresponding to 
the exact soliton complex for a = 3/4. To investigate the possible existence 
of other forms of soliton complexes, let us insert the soliton anzats (31) with 

uo{z) ~ 8arctan(exp(2;//)), (44) 

into the potential expression Eq. (43). Then the linear equation (42) takes the 
form 

After introducing the new coordinate y = z/l Eq. (45) can be rewritten as 

l-ar'^^, + n-r'-^m)-o^ (46) 

Solutions of Eq. (46) give us the next approximation for u^z- 

Recalling asymptotics of the eigenvalue problem and Eq. (40), we have to 
impose the following condition on parameters a and I: 

ar^ = 1 + h - r^ = e, (47) 

which must hold for any soliton complex. We have introduced the parameter 
e which takes, evidently, a discrete set of values e„ . Thus finally we find the 
equation: 
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Eigenvalues and eigenfunctions of this type of the equation are well- known [47] . 
In particular, the equation for parameter e„ reads: 

l = i[(l + 86-1)1/2 _2n-lp. (49) 

As a matter of fact, the problem is now reduced to the determination of all 
values of the potential depth in Eq. (48) for which the discrete level equals 
unity. From Eq. (49) it follows 

e„ = 2/[(n+l)(n + 2)], (50) 

where n = 1,2, 3.... Then one obtains for Z„ and q;„ 



where cr„ equals 



and 



C2 = /i + cr„, (51) 

n{n + 3) 



(n+l)(n + 2)' 



(52) 



a„=a^(l + — )-2. (53) 

Here a° are eigenvalues of the problem for the case h = 0, i.e. for the dispersive 
sine- Gordon problem: 

^ 2(n+l)(n + 2) 

[n(n + 3)]2 ^""^^ 

Ten first eigenvalues a^ are shown as solid circles in Fig. 4 (to show more vividly 
the dependence on n we connect the points by the solid line) . The infinite series 
of discrete values of a° rapidly diminishes with increasing n. Only odd values of 
n seem to be valid for the problem under the consideration. However, as shown 
in Fig. 4, the even values of n also reproduce well enough the exact eigenvalues 
(solid squares) which we find by numerical integration of Eq. (16) (see section 
6). This is explained by the fact that the next iteration step leads to splitting 
of levels with n > 1 into two close levels with even and odd eigenfunctions. The 
level splitting is smaller for higher n, and the analytical dependence a^ serves 
as a good approximation of the eigenvalues in the dispersive sine-Gordon case. 
Eigenfunctions of Eq. (48) also describe well the changing of the form of 
47r-soliton when it appears in the excited state. Functions tpn, corresponding to 
Uzz, vanish exponentially at z = ±oo and exhibit an oscillating behavior at the 
coordinate origin. The oscillation domain increases for higher n-values. This 
oscillation can be interpreted as the radiation locked between two 27r-solitons. 
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Thus the sohton complex can be reahzed by the infinite series of configurations 
to be referred to, naturally, as "excited" states. 

Figure 5 presents a„ of five levels as functions of the parameter h. We see 
that values of a„ tend quickly to zero at high fields, reflecting correctly the qual- 
itative tendency in the /i-dependent behavior of exact eigenvalues. Quantitative 
comparison will be done in section 6. 



5. EXACT SOLITON COMPLEX SOLUTIONS IN DOUBLE 
QUADRIC MODEL 

Explicit expressions for 47r-soliton complex and its "excited" states can be 
found in the framework of the double piecewise-linear model Eqs. (13) and (14) 
[48] . In this case stationary moving complexes can be constructed in analytical 
form , and hence, the corresponding eigenvalue problem for the parameter a„ 
can be exactly solved. 

In this section we exhibit principal results of the considerations. For the 
sake of simplicity we present main formulas in the limit case ft, = in reason of 
their clarity. 

As in the dDSG case, after introducing the coordinate z in the moving 
reference frame, we derive the equation: 

Uzz + auzzzz - .f{u) = 0, (55) 

where f{u) is given by the expression Eq. (14). We are interested in odd 
solutions of Eq. (55), u{—z) = —u{z), with limit conditions m(±oo) = ±27r. To 
construct the 47r-soliton it is sufficient to find the solution in the external region 
(z > zo) and the general odd solution in the internal region (|z| < uq). They 
look like the foUowings: 

Ue{z) — 2Tr — A C'Xp{—K2z), (56) 

and 

Ui{z) — B sm(kiz) + Csinh(Kiz), (57) 

In the limit case h = the exponent K2 coincides with the parameter ki of the 
solution (57): 



«i = «2 = AC = ±\/[(i + a)i/2 _!]/«. (58) 



The parameter fci equals to: 



A:i = fc=\/[(i+a)V2 + i]/a. (59) 
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Using the conditions of continuity of the function u and its first, second and 
third derivatives in the point zq, where w — mq, we arrive at the following set of 
equations: 

i?sin(fczo) + C'sinh(KZo) = 27r — j4exp(— kzq) ~ uq, (60) 

fci?cos(fczo) + KCcosh(KZo) = AKexp(— KZo), (61) 

— k _Bsin(fczo) + K Csinh(KZo) = — Ak exp(— kzq), (62) 

— fc'^i?cos(fczo) + tC' C cos}i{kzq) — Ak'^cxp(— KZo). (63) 

From the system Eqs. (60-63) one can find the eigenvalues of the parameter a 
and the coefhcients A, B and C. 

It is easy to see from the Eqs. (61) and (63) that cos(fczo) ~ 0, and hence 
zq is given by 

TT 1 

za = -{n- -) = fi,Jk. (64) 

Introducing the parameter X{a) = K{a)/k{a) we obtain the following equation: 

A = exp(-A^„), (65) 

which determines the eigenvalues A„ and a„. It is clear that this equation has 
a solution for every n. Moreover it is reduced by the substitution A„ = A/i„ to 
the Lambert's equation: 

A„exp(A„) = ^„. (66) 

Its solution is the Lambert's function A„(a) = W{nn)- By solving the last 
equation with respect to a we find that there is a infinite series of the a„- 
values for which we can construct the 47r-soliton. Coefficients of the solution 
are expressed through fc(Q:„) and K{an) as follows: 

A = 7r-, (67) 

B = i-ir-',S^, (68) 



and 



^ 2TrkK , , 
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Ten eigenvalues a„ are shown as the open circles in Fig. 4. Note that 
parameters a„ diminish with increasing n in the like manner as in the analytical 
results Eq. (54). 

In the general case ft. 7^ the equation for a„(/i) is much more complicated 
[48] than Eq. (65). The results of its solving are presented in Fig. 6. One can 
see that a„(ft.) are quickly decaying functions of h. 

Now we discuss properties of the cigenfunctions. Evidently, the soliton com- 
plex states can be classified by the integer values of the parameter n. In the 
"excited" complexes the value of n shows the number of nodes of the oscillating 
part of the solution. The latter corresponds to the radiation locked between two 
composite solitons. 

The first five cigenfunctions are shown in Fig. 7, where they are numbered 
from the left to the right by n = 1,2,3,4,5, respectively. For convenience of 
observation we have shifted the centers of complexes in space. With increasing 
n the separation L between the 27r-solitons slowly grows. For large n, L ^ 2zq ~ 
21n(7rn/ln(n)). At the same time the amplitude of the oscillations decreases, 
and in the limit n ^ 00 the soliton complex is approximated by two well- 
separated solitons and a linear standing wave between them. 

It should also be noted that Uz>Q for all functions. More precisely, at the 
soliton center Zc the first derivative reaches its local maximum value for an odd 
n and Mz = for the even n. The latter is possible if the expansion of u{z) 
at the point Zc begins with a term of the order (z — ZcY . Indeed, it is easy to 
be convinced that the condition kB + nC = is fulfilled for the solution (56), 
(57) with the even n. The same behavior of the cigenfunctions was found while 
integrating the dispersive sine-Gordon equation [11]. For the general case h ^ Q 
the condition Uz{zc) — does not hold in both the piecewise- linear and dDSG 
models, and the first derivative becomes positive at the soliton center for all n. 

Thus the simple piecewise-linear model exhibits main peculiarities of soliton 
complex structure and its stationary dynamics, which turn out to be universal 
for the dispersive nonlinear models. 

6. SOLITON COMPLEXES IN THE DISPERSIVE DOUBLE 
SINE-GORDON EQUATION 

In this section we present results of the numerical intergation of the dDSG 
equation in the case of the stationary complex motion. We start with Eq. (16) 
and seek its 47r— soliton solution, i.e. impose the limit conditions: 

m(— 00) = 0, u(oo) = 47r, ?iz(±oo) — ^^^(±00) = Uzzzi'^oo) — 0. (70) 

However, at first, we propose one more effective approach to solving the 
equation. It allows us to simplify the integration procedure and formulate some 
strict assertions about a possibility of different forms of the soliton complexes. 
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Eq. (16) describes the effective particle dynamics in the four-dinicnsional 
space {u,Uz,Uzz,Uzzz}- It has the first integral: 

/ = a[{ul)zz - 3mL] +ul- 2(1 - cos{u)) - 8h{l - cos(|)). (71) 

By virtue of the hmit conditions, Eq.(70), / = in the case of the soliton 
complex solution. In general, the phase space is determined by four variables. 
However if we found the solution u{z), we could obtain the trajectory in the 
two-dimensional phase space, {u,Uz} by expressing Uz as a function of u. 

Let us introduce the function u1 — F{u) or Uz = \F ■ By the differentiation 
of the definition we get that Uzz — \Fu. After substitution of these expres- 
sions to Eq. (71) we find that the function F obeys the following second-order 
differential equation: 

1 u 

F + a{FF^u - -^Fl) = 2(1 - cos(w)) + 8/i(l - cos(-)). (72) 

From the conditions of Eq. (70) it follows that F{u) must be periodic with the 
period At: and i^„ = for m = 0, 47r, .... It is easy to find the function expansion 
at u = 0: 

F^K^u^ + 0{u'^) + ..., (73) 

where k is given by Eq. (40). 

To this end we are able to perform the numerical integration of the equation 
(72). This is achieved by use of a fourth-order Runge-Kutta method. To study 
the dependence of the soliton complex form on the parameter h we begin the 
integration, at first, for /i = (the dSG case). The corresponding results are 
presented in Figs. 4, 9, 10(a), 10(b). In this case some results can be also 
verified and compared with those found in the paper [11]. 

First of all, it is easy to be convinced that there is no 27r-kink solution in Eqs. 
(72) and (16). Next, there are only discrete set of q;„ for which the 47r-soliton 
complexes exist. We show the first ten a„ as the solid squares in Fig. 4 and 
see that they are in a good quantitative agreement with the analytical results 
(54) (Fig. 4, solid circles). In Figs. 10(a) and 10(b) the phase portraits are 
drawn in the plane [u^Uz)- The main soliton complex trajectory is realized for 
ax — J and has only one maximum. Other phase portraits differ by the number 
of maxima. Call this n. So it is convenient to classify these states by the integer 
n using the direct analogy with the results of phenomenological approach and 
the piecewise-linear model. 

In particular, as in the case of the piecewise-linear model all the odd states 
(Fig. 10(a)) have the nonzero first derivative at the soliton center where u = 
27r. The functions corresponding to even n (Fig. 10(b)) demonstrate another 
behavior. Their first derivatives at the point u = 2it behave as u^ ~ (u — 27r)^'^. 
The reason of the specific dependence is the same as that in the piecewise-linear 
model, namely, it turns out that u{z) — 27r ^ (z — z^)"^ and hence Uz ^ {z — Zc)^- 
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As a result Uz — in the point u — 2tt for states with the even n. However as 
we show below, the inclusion of the magnetic field, i.e. taking into account the 
term with nonzero h, removes this degeneracy. 

Now the functions u„(z) for the complex and its "excited" states can be 
found by numerical integration of the equation u^ ~ \/ F{u). First five soliton 
complex states are presented in Fig. 9. They virtually have the same shapes as 
solutions of the piecewise- linear model (c.f. Fig. 7). 

When h ^ the new mechanism of the soliton attraction begins to work. The 
magnetic field draws together the composite solitons. Results of the numerical 
solution of the Eq. (72) for this case are shown in Figs. 11, 12(a) and 12(b). 
One can see that the phase portraits of all the " excited" states become similar 
with increasing n even at small h. This means that in this case the form of the 
"excited" complex approaches the wobbler solution. It is clearly seen from Figs. 
12(a) and 12(b), where solutions for h = 0.1 are presented. It also confirmed by 
the analysis, performed below, of the dynamical characteristics of the complex 
and its "excited" states. 

Now, we discuss the dependence of q;„ on the parameter h. For the n — 1 
there are exact soliton-complex solution of the dDSG equation and the analytical 
expression for ai (h) is given by: 

a^{h) = l{l + ^hr^. (74) 

For the "excited" complex states the dependences an{h) are found numerically, 
and the first five of them are presented in Fig. 8. They are qualitatively similar 
to those of the picccwise-linear model (Fig. 6) and of analytical results (Eqs. 
(52)-(54) and Fig. 5). Quantitative comparison reveals that with increasing 
h numerical eigenvalues a„{h) vanish more rapidly than the analytical depen- 
dences. However we believe that the decaying of the functions at high h is 
proportional to 0(1 /h^). Wc have fitted the data for a„(/i) by the expressions 
of the form Eq. (53) 

"n-""(l + 7^)-', (75) 

where Z?„ is the only fitting parameter. As a result we have found a good 
coincidence between the data and the analytical approximation. For example, 
in the case n = 2, the value D2 ~ 2/5 provides a deviation less than some 
percents on the interval < h < 2. In general, £*„ decreases quickly with the 
growth of the number n. 

At high h the soliton complex and wobbler shapes have to degenerate to the 
kink of the sine-Gordon equation for the variable (j) = u/2. Indeed, introducing 
this new variable one can exactly rewrite Eq. (16) as the following form of the 
stationary dispersive double sine-Gordon equation: 

(/)^^ + q:(/'55{{ — sin(0) — hsiii((f)) cos(0) — 0, (76) 
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where the following notations are introduced: ^ = yhz^ h — l/h, and a = ah. 
When (5 = 0, Eq. (76) has the wobbler solution which is reduced to 27r-kink in 
the limit h ^ 0. The soliton complex and its "excited" states exist at /i 7^ and 
Of ^ 0. As it follows from Eq. (75), at high h the linear dependence between 
parameters a and h is suggested, i.e. 



a. 



„0 n2 



DlK (77) 

In particular, in the case of the exact solution we find from Eq. (74) that 

&i{h)^^{l + lh)-^ (78) 

and, at small /i, in fact, ai « h/3. The eigenfunction 0i(C) would be considered 
as the direct continuation of 27r-kink solution of the SG equation for the case 
h^O. 

The eigenfunctions '/iti(C) for n > 2 behave like the wobbler solution. To 
understand the proximity of the " excited" states and the wobbler we apply the 
iteration procedure (see section 4.B) to Eq.(76), starting with the wobbler-like 
ansatz (31): 

0(0 = 2 arctan{exp((j^ - R)} + 2 arctan{exp((7^ + R)}, (79) 

where now q and R are the constant parameters of the solution. The ansatz 
seems to be advantageous over the previous one, Eq. (44), because it includes 
two parameters. However, one must keep in mind that already after a first 
iteration step all parameters get definite values, and the eigenfunction form will 
be corrected after every step. Here we use the appropriate choice of the ansatz 
to obtain the analytical estimation for the soliton separation R in the complex 
at high field h. 

Omitting details of the calculation we derive finally the following equation 
for the determination of an{h): 



dte cosh2(e)-sinh2(i?)^ 



6iq'-T-:^ + l + h-q^ Vl r^T^]'^« = ". (80) 



From this equation, in the case of small h, we find that the quantity Uq = 
2{h — sinh^(_R))/ag^ has to be equal to {n+ l)(n-|-2). It is easy to see that this 
relation is in a qualitative agreement with the assumption about the behavior 
of din (see Eq. (77)). Then we can use Eq. (77) and put q w 1 to obtain the 
following relation for the determination of the parameter R: 

smh^{R)Kh{l-Dl). (81) 

As it follows from Eq. (81), when n increases, the fast decay of values £)„ causes 
the parameter R rapidly approaches the wobbler value R^ (see Eq.(15)), even 
at small h. 
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Thus, the smahness of a„(ft.) and its rapidly- vanishing dependence on h and 
n, lead to the degeneracy of the complex "excited" states to the wobbler-like 
solution. 

This fact can be demonstrated perfectly, if one calculates the energies of the 
soliton complex and its "excited" states as functions of the parameter h (see 
Fig. 13). We have normalized these quantities by the energy of the wobbler. 
As clearly seen from Fig. 13, energies of the "excited" states approach the 
wobbler energy very quickly with increasing h, while the exact solution remains 
well-separated from the wobbler up to high h. 

Another form of the dDSG equation (76) can be useful to analyze the hi- 
erarchies of the bifurcation values of parameters am{h) for other multisoliton 
complex solutions. One can note that there is an infinite series of solutions of 
Eq. (76) for h — Q and a = a^^ corresponding to the 47r-complex for the vari- 
able (/). When ft, 7^ these solutions are modified, giving new branches am{h). 
It is evident that the eigenvalues correspond to the 87r-soliton complex and its 
"excited" states for the variable u{z) obeying Eq. (16). It is clear that at high 
h the behavior of eigenvalues am(h) is approximated by a^/h. This is a quali- 
tatively different behavior than that in the case of the 47r-complex. At small h 
the parameters am{h) reach their constant bifurcation values 0:^(0), which can 
be found from the numerical solution of Eqs. (16) and (71). 

It is important to note that a concrete physical system is characterized by 
the definite value of the dispersive parameter /3. Then for the given /3 and h 
several soliton complex states can exist simultaneously. 

In the case of the IdDSG equation at the fixed (3 there exists a discrete set 
of velocity values Vn{h) corresponding to the stationary radiationless motion of 
the complex: 



Vnih) - y^l - (/?/«„ (ft)) 1/2, (82) 

where a„{h) are eigenvalues of the nonlinear spectral problem (16). As an 
example, the dependences of 14, (ft) for the soliton complex and its "excited" 
states at /3 = 1/12 arc pictured in Fig. 14. Every branch has a finite range of 
the velocity changing from the maximum values V^ = ^/l — {/3/a^)^/-^ to zero. 
It is evident that critical fields ftn corresponding to static solutions are found 
from the equation Q!„(ft„) = (3. The number of possible complex states is finite 
because the states with Q:„(ft„) < /3 turn out to be forbidden. 

The regularized 2dDSG equation (10), at any /3, has a complete infinite series 
of soliton complex states with the following allowed velocities: 



In the case of n = 1 we have the explicit expression for ai{h) (Eq. (74)), and 
after its substitution to Eq. (83) we come back to Eq. (20). It is clear that the 
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maximum velocity value is realized at h — 0, when a — 3/4; it coincides with 
Vr{f3) of Eq. (7). At high h and hence small an{h) the velocity dependences 
are simphfied to Vn{h) ~ A/a„(/i)//3. 

In conclusion we present the energies of the complex and two " excited" states 
as functions of h for a prescribed /3 = 1/12 (Fig. 15). As shown above (see Eq. 
(21)), El is independent of h. Energies of the "excited" states turn out to be 
less than Ei; therefore this terminology becomes inadequate when (3 is fixed. It 
seems also that the energy arguments about the complex stability do not work in 
this case. However, it should be noted that the velocities of all the complex states 
are different, and the exact complex solution possesses the highest velocity. A 
stability criterion has to be formulated so that the energy comparison is to be 
performed under the condition of the conservation of the another integral of 
motion, the momentum. Since for a fixed /3, both first integrals are functions of 
the parameter h alone, then for a given h, the momenta are different, so that 
such a simple comparison becomes impossible. The understanding of the above 
studied mechanisms of the formation of the topological soliton complexes allows 
us to believe in their stability in the dispersive conservative systems. 

7. SUMMARY AND DISCUSSION 



In this section we discuss some peculiarities of soliton-complex dynamics in 
strongly dispersive models and formulate final general conclusions. 

Dispersive nonlinear equations with fourth-order spatial derivatives, exam- 
ples of which have been studied above, are used usually as a first approximation 
in the description of discrete systems. The present investigation shows that the 
higher dispersive terms in the nonlinear wave equations effect much more essen- 
tially than the small perturbations. The dispersion causes the strong dissipation 
of energy of the moving 27r-soliton, but also makes it possible the creation of 
the bound soliton complexes consisting of two or more 27r-solitons which can 
move radiationlessly. 

The formation of soliton complexes turns out to be the universal phenomenon 
in nonlinear strongly dispersive media, which are of both theoretical and prac- 
tical interest. 

From the theoretical point of view it is very interesting that the dispersion 
produces a discrete spectrum for a nonlinear eigenvalue problem for purely soli- 
tonic solutions, and its influence is not reduced simply to small changes of the 
27r-kink shape as it was believed before. Analytical and numerical considera- 
tions confirm the absence of the steady 27r-kink solution in the above-studied 
dispersive systems with the fourth and higher order spatial derivatives. The 
earlier papers (see the review [49] for references) declared the existence of such 
solutions. Indeed, application of the perturbation theory, using the smallness of 
a, leads to this erroneous conclusion. It seems to be possible to construct for- 
mal asymptotic series for the 27r-kink solutions of Eqs. (16) and (72). However, 



23 



summation of all terms of the series results in disappearance of the periodicity 
of the function F{u) with period 27r, and only solutions with periods AtTjGtt... 
survive. 

At the same time these peculiarities of the dispersive effects require some 
caution when one exploits the variational description of the bound solitons or 
breather states in dispersive or discrete systems [50] . 

The numerical finding of the discrete set of periodic eigcnfunctions in the Eq. 
(72) could point to the possibility of a complete integrability of the equation, 
at least, in some special cases. However this question is still open. 

The existence of soliton complexes consisting of three or more solitons is not 
discussed in detail in the given paper. The above analysis of the mechanisms 
and conditions of the creation of the complexes is applicable to the case of the 
multisoliton bunching. The presence of exact solutions in Eqs. (24) and (27) 
may be used as a basis for further studies in this direction. One is easily con- 
vinced that the exact solutions found are not exotic. Specific relations between 
constants (3 and 7 of Eqs. (23) and (26) are required only as conditions of the 
existence of the solution in a very simple analytical form. If we slightly alter 
the value of the parameter 7, the solution still exists, just as before. Numerical 
integration of the dSG equation in the stationary case [11] revealed such a so- 
lution for 7 = 0. Hence the multicomplex solutions are stable with respect to 
changing the dispersive parameters. 

The important question of the stability of the soliton complexes in the frame- 
work of the theory of partial differential equations has remained beyond the 
scope of the paper. Work in this direction is in progress. However, previous 
studies of the problem in the discrete and nonlocal sine-Gordon models suggest 
a positive answer to this question. 

From the experimental point of view the soliton complexes may be very at- 
tractive for the use in energy and information transfer processes. In particular, 
the discrete arrays of Josephson junctions can be described by one-dimensional 
dispersive sine-Gordon models [30,31,50]. The fiuxons, the 27r-kinks, in these 
systems could form the soliton complexes with specific properties discussed 
above. Therefore, Josephson junction arrays may be considered as suitable 
candidates for real experiments, in which the soliton complexes would mani- 
fest themselves [15]. Other examples of appropriate physical systems are the 
low-dimensional ferro- and antiferromagnets. More definitely, it may be an 
one-dimensional biaxial ferromagnet with strong easy-plane anisotropy and a 
magnetic field applied along the plane. This system is described by the dou- 
ble sine-Gordon equation, and taking into account the discreteness effects duly 
leads to the dispersive model (9). Two-dimensional antiferromagnets with the 
weak interplane exchange can be also treated in the framework of the models 
analyzed in previous sections. The nonlinear dynamics of dislocations is another 
potential field of application of the obtained results [1,2]. Thus there are many 
experimental possibilities for the observation of the occurrence of stationary 
soliton-complcx dynamics in dispersive media. 
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In conclusion we note that the present investigation shows that purely soli- 
ton complexes can be realized in dispersive media which are modulationally 
unstable from the very beginning, as well as in systems where the 27r-kink mo- 
tion produces the medium modulation due to the radiation accompanying the 
soliton movement. In both situations nonstationary radiative dynamics of two 
solitons can result in the final formation of the purely solitonic complex with 
cancelling the oscillations on soliton wakes, what can be considered as a specific 
interference efi^ect. However, discussion of the nonstationary soliton complex 
dynamics as well as the problem of their dynamics in dissipative systems is the 
subject matter of following publications. 
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APPENDIX 

A. Effective Lagrangian and Hamiltonian 

After inserting the anzatz (31) into the Lagrangian (11) and the Hamiltonian 
(12) we perform all needed integrations using the following values of intergals: 



— oo 

h{R) = 7 dy{r^^2Mm£2Sh^)2 ^ 1+ . gi? 
^ ^ J^ yncosh2(fl)+sinh2(y)/ smh(2i?) 

^.(^) = 7jyi,^fmm^,r - coth^(i?)(i - ^41^) 

— OO 



-oo 



As a result expressions for the kinetic and potential energies of the two-soliton 
system can be written as 
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T = ^{llt^[2R^ + (i?2 + \)h{R)] + {lRtt^nh{R)fl2{R) + 2RRtllt 

+Xt^h{R)}, 



(84) 



and 



U = f{Ii{R) + PhiR) + 2hPloiR) 
^f[l-tanh^(i.)(l-/.(i?)(l + -^))]}. («^) 

Then the effective Lagrangian and Haniiltonian are expressed as usual 

L^ff^T^U, H,ff^T + U. (86) 

The Lagrangian function with accuracy to order 0{R^) can be written 

L = U + R^[Ti-Ui{q,Xl)l (87) 

where the zero-order Lagrangian is the foUowing one: 



¥4?— -<^")'^-4 



-{-li-l + Xt-{- + Kr + ^}, (88) 



and Ti contains all terms proportional to the time derivatives of function I and 
Ui is given by: 

t/i(g,X,2) = i^{X,2-l + ^ + lf}. (89) 

From Eq. (87) we see that one of the Lagrange equation takes the form dL/dR = 
0. It is transformed into RUi{q, Xf) = for stationary states. 



B. Intergable sine- Gordon limit. Exact two-soliton solution 



If we put h = 0, /3 ^ 0, I ~ J^^iv) = Vl — v"^ = const, and Xt = 0, then 
we reduce the problem to the consideration of the sine-Gordon system in the 
reference frame moving with mass center. In this case the Hamiltonian i?e// is 
simplified to Hsg- 

HsG = ^^^^ {(1 + {t'<^nh{R)Rtf)h{R)i-^ + h{R)}. (90) 

It is easy to be convinced that the function R{t) 

R{t) = ln[icosh(7i;i) + W (-cosh(7wt))2 - 1] (91) 
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is the exact solution of the effective Hamiltonian equations. The substitution 
of R{t) into the ansatz (3f ) and Hsg provides the exact two-sohton solution of 
the intcgrable SG equation and corresponding value for the energy Ea — I67. 



C. Double sine-Gordon limit. Small wobbler oscillations 

Suppose again that (3 = but let the parameter I be the function of R of the 
form 1{R) ~ tanh(_R). Then after substitution of 1{R) into the Lagrangian we 
find for the double sine-Gordon system: 

Lose = SRt^ tanh(i?)G(i?) - coth(i?) - hR. (92) 

The last two terms in Eq. (92) correspond to the potential energy, and the 



function G{R) is introduced as 



TT^ 



G{R) = h{R) + —-^^^^[2R' + h{R){R' + -)]. (93) 

The stationary solution is found from the condition of minimum of potential 
energy. One determines the equilibrium value of the parameter R from the 
equation sinh(i?^) = 1/v/i. Naturally, the ansatz coincides with the exact 
wobbler solution (15). 

Small oscillations near the equlibrium position r{t) = R(t)—Rw are described 
by a linear equation which is found from the Lagrangian by the usual way. 

nt + n\ = 0. (94) 

The frequency is given as 

n^(h) = 2hcoth^{R^{h))/G(R^(h)). (95) 

This field dependence for the frequency of the internal oscillations is seemed 
to be the best analytical approximation of this function for the moment (see 
[51] and references therein for the comparison). 
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FIGURES CAPTIONS 

Fig. 1. Velocities of soliton complexes as functions of the discreteness pa- 
rameter d. Solid (dash) line corresponds to the IdSG (2dSG) equation. 

Fig. 2. The potential shapes of the dispersive double sine-Gordon and 
piecewise-linear models. 

Fig. 3 (a)-(d). The potential energy of the dispersive double sine-Gordon 
system in the variational approach as the function of the effective length and 
the separation between solitons. Figs, (a) and (b) correspond to (3 = 1/5, and 
(c) and (d) to /? = 3/4, respectively. The parameter ft, = for the cases (a) and 
(c), and h = 0.1 for the cases (b) and (d). 

Fig. 4. Eigenvalues «„ {n = 1..10) corresponding to the soliton complexes. 
Solid circles denote results of the phenomenological approach to the stationary 
dispersive sine-Gordon equation. Open circles correspond to the piecewise-linear 
model. Solid squares are the eigenvalues obtained by the numerical integration 
of the dSG equation. 

Fig. 5. Analytical dependences a„(ft) (n — 1..5) found after one step of the 
phenomenological iteration procedure for the dDSG equation. 

Fig. 6. Dependences an{h) (n = 1..5) obtained for the dispersive piecewise- 
linear model. 

Fig. 7. Five first soliton complex eigenfunctions constructed explicitly in 
the piecewise-linear model (the case h = 0). 

Fig. 8. Numerical results for the eigenvalues a„(/i) (n — 1..5) in the dDSG 
equation. 

Fig. 9. The soliton complex and its first four "excited" states found numer- 
ically for the dSG stationary equation. 

Fig. 10 (a),(b). "Phase portraits" of soliton complex and the odd (a) and 
even (b) "excited" states in the dSG equation. 

Fig. 11. The soliton complex and its three "excited" states found numeri- 
cally for the dDSG equation (h = 0.1). 
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Fig. 12 (a),(b). "Phase portraits" of soliton complex and the odd (a) and 
even (b) "excited" states in the dDSG equation (ft, = 0.1). 

Fig. 13. Energy dependences of the soliton complex and the " excited" states 
on the parameter h. The energies are normalized by the wobbler energy E^ih). 

Fig. 14. Velocities of the soliton complex and the "excited" states as func- 
tions of the parameter h. The dispersive parameter /3 equals 1/12. 

Fig. 15. Energies of the soliton complex and two "excited" states as func- 
tions of the parameter h. The lines finish at critical values /i„. The value 
hi — 4/3 is beyound the figure domain. The parameter (3 is fixed as 1/12. 
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